o ■ Effect of symmetry of the electron states of HTSC on the 

^ current-voltage characteristics of SIS junctions. 

u ; 

<d : P.I.Arseev., N.K. Fedorov, S.O. Loiko 



> 
in 



P.N. Lebedev Physical Institute, Russian Academy of Science. 
Leninskii pr. 53, Moscow, 119991 Russia. 



cn 

o 

o _ 

^ I Abstract 

^ ■ The current- voltage IV characteristics of SIS junctions are calculated in the frame- 

, work of a multiband model with an anisotropic effective order parameter of HTSC. The 

"j^ I results of calculations show that the shape of IV characteristic and the density of elec- 

tron states changes significantly depending on the parameters of the model. A theoretical 
, . explanation is proposed for the experimentally observed s-like behavior of the IV charac- 

^ I teristics of SIN and SIS junctions with BSCCO-type superconductors. The dependence 

O ' of the superconducting peaks asymmetry on the mutual arrangement of the bands is an- 

^ ■ alyzed. The difference between the obtained results and the results of single-band models 

with the s and d symmetries of the order parameter is discussed. 

The lack of a generally accepted explanation of the pairing mechanism in high-temperature 
CN '. superconductors (HTSC) makes the interpretation of numerous experimental data (including 

in 



the results of tunnel experiments) difficult. The complexity of the crystal structure of HTSC 
O compounds also complicates the situation. In most cases, it is difficult to explain peculiarities 
2 ■ in the tunnel characteristics within the framework of the standard Bardeen-Cooper-Schrieffer 
(BCS) model. Such peculiarities include the variety of the subgap structure ranging from 
^ , a linear d-type structure to a nearly planar structure of the s type P, the asymmetry of 
G [ superconducting peaks at the edge of the gap [3 0], and nontrivial behavior of the tunnel 

^ ' density of states of SIN junctions and of the current- voltage characteristics of SIS junctions 

Ch dl 

O outside the gap region. In some experiments [H El , strongly suppressed —— characteristics with 
O : dV 

very narrow peaks typical of s-type superconductors were obtained in the subgap region, which 

• I— I contradict at first glance the available data on strong anisotropy of the order parameter in the 

Cu02 plane f6J. 

\ A large number of theoretical models have been proposed to explain the variety of experi- 
mental data. These include the models based on various features of the band structure (taking 
into account, for example, the closeness of the Van Hove singularities to the chemical potential 
[7, 8 ), on the choice of the order parameter symmetry jH], and on the properties of the tunnel 
matrix element jH] in the single-band BCS scheme. Even the slave-boson approach JU] was 
used for explanation the asymmetry of the superconducting peaks and the inhomogeneity of 
the density of electron states in the normal state. A series of publications IT^ IT!^ ITl] is 
devoted to the role of inelastic scattering, including the inelastic scattering by antiferromag- 
netic spin fluctuations. The results of theoretical investigations show that the inclusion of Van 
Hove singularities and the use of the (s -|- ci) symmetry of the order parameter are obviously 
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essential for a correct explanation of the most experimentally observed properties of the tunnel 
conductivity of HTSC For this reason, the experimental investigation of the properties of the 
electron spectrum of high-temperature superconductors in the normal state is very important. 

In our opinion, it is the analysis of the crystal structure that makes it possible to explain 
a number of results of tunnel measurements in high-temperature superconductors. This paper 
aims at explanation the features of tunnel characteristics using the model with an anisotropic 
effective order parameter, which is based only on the properties of the electron spectrum of 
HTSC jin]. The spectrum of the considered model approximately corresponds to the band 
structure of compounds of the BiSrCaCuO (BSCCO) type. An important feature of this ap- 
proach is that the initial electron-electron interaction leading to pairing is treated as isotropic 
in the Cu02 plane. The anisotropy of the order parameter and of the excitation spectrum is 
determined by the symmetry properties of the crystal lattice. According to our calculations, 
many features of the tunnel characteristics typical of HTSC compounds can be explained by 
using a rather universal approach which will be described below. 

Let us consider the effect of symmetry of the initial bands of the superconductor on the 
current-voltage characteristics of S'JS'-junctions. The electron system of the Cu02 plane is 
described by model Hamiltonian of the form: 

i,cr i,cr i.cr 

+ E it%pA^P-(j^ ^) + ^-c-) + E it'z'-py4,aPyij, ^) + h.c.) + 

+ E (titpJtPxij, (t) + h.c.) + J2 {t'lp/tPyij, + h.c.) + 

+ Uz2j2cf^Ciict^Ci^ + UdJ2dtidiidt^di'i , (1) 



Here cf^^ and df^^ are creation operators for electrons with spin a on the c?^2- and d^^-yi- orbitals 
of the i-th. site of copper; a) and Py{i^ a) are the creation operators for electrons with spin 
a on the Px and Py orbitals of the i-th. site of oxygen (Fig.l). e^, and £^2 are the energies 
of the p levels of oxygen and of the d^^-yi and d^^ levels of copper, respectively, measured 

from the chemical potential (/i = 0). t^J_p^ = ^^S-py = ^z^~p and t^^-p:, — ~^d-py = ^d-p are 
the matrix elements of one-particle transitions between the dz2 and d^-i-y-i orbitals of copper 
and the p orbitals of oxygen. For the sake of simplicity, we assume that superconductivity is 
caused by the isotropic attraction f/22 < of electrons on the dz'^ orbital. It should be noted 
that a more complicated consideration of the superconducting interaction would not lead to 
a qualitative difference from the case investigated here. We also assume the presence of the 
isotropic effective electron-electron interaction Ud on the dx2_y2 orbital of copper. 

Introducing new operators a[f]. (/3=1,2,3) in the k representation according to the formulas 

^ _ I r^(2)^(2) , ^(3)^(3) 



p.(k,a) = C7£)a«+C(f)aS + C£)a£ 



(2) 



PyiKa) = Cl^Jail + Cl^Ja^l + Cl'Ja^l, 
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Figure 1: Schematic representation of atomic orbitals in the CUO2 plane and transitions 
between these orbitals included in Hamiltonian (1). The inset shows schematically the arrange- 
ment of energy bands and levels corresponding to Hamiltonian (4). 



we diagonalize the Hamiltonian part describing one-particle transition between the dx^-y^ or- 
bitals of copper and the p orbitals of oxygen. The coefficients and Cp^i^^) in formulas (2) 
have the form: 



^(1) _ n /^(") _ 



cy> = 0, cr = 



^{ep-s^y + itxr + ity 



'^^''^ ^/{Q' + {tyf ''^'^ ^{ep-e^Y + itxY + ityf 

where o; = 2,3 and tx(y) = t^-p sin . Similarly, the matrix element of a transition between 
the (i;j2 orbitals of copper and the p orbitals of oxygen in the k representation can be written as 
^x{y) — iz2_pSin Since the matrix elements and t^^i_p depend on the indices of different 

sites, the quantities tx{y) and t^^y) must be functions of k. Operators , a^^^ and a^^^ are 
the creation operators for electrons with spin a and quasimomentum k in the bands formed by 
hybridized p orbitals of oxygen and dx^^y2 orbitals of copper with the dispersion relation 



s, = ep , e2(3)(k) = T y{sp - + A{tl{k) + ^^(k)) . (3) 

As a result of transformation (2), the system Hamiltonian in the mean-field approximation 
takes the form: 

+ E W^XA'l + h.c.) 

k,(T,/3 



k k,/3,7 



where 



W,{k) = -^-^^^Iht , W^{k)= , ^^^^ a = 2, 3 (5) 

and new order parameters 

A^,(k) = Cf (k)ci^Vk) Ad , 7 = 1, 2, 3 , (6) 

have been introduced, in which 

= -^t^dE(^k,TC^-ka) > = -^U,2 ^(Ck,tC_k,i) • (7) 

The arrangement of the energy bands corresponding to Hamihonian (4) is shown in the inset 
to Fig.l. In the framework of a reahstic description of BSCCO-type compounds it is 
assumed that the chemical potential /i and the z"^ level of copper locate closely to the center 
of the upper band £:3(k). The remaining parameters of the model are assumed to satisfy the 
relations Ed < Sp, Sp- td-p, t^2_p < |ed|, |Ad| < A^2 < \ed\, \ep\. 

Since only the energy range near the chemical potential (on the order of several A22) plays an 
important role here, we can disregard the two lower bands £1 and £2(k) taking into account the 
relation between parameters and restrict our analysis to the upper band e3(k) only. In addition 
to one-particle transitions included in Hamiltonian (1), we can take into account the additional 
hybridization of the dz-2. orbital with the atoms surrounding the CUO2 plane, including BiO 
and SrO complexes of the BSSCO-type compounds as well as direct transition between the 
dz2 orbitals of different copper atoms. In this connection, we assume that the initial band is 
characterized by the dispersion relation e^^ik) = £^2 +tz2(cos{kx) +cos{ky)) with a width much 
smaller than for e2{k.) and £3(k) bands. Thus, in the chosen approximation, the initial problem 
can be reduced to a two-band model in which the initial (k) and the e3(k) bands with the 
one-particle hybridization 14^3 (k) and the order parameters Az2 and A33(k) are considered. 

We introduce the following time Green's functions: 



GAk;t,t') = -i(T c^,,it)cl„it')) , w^.t(k;t,t') = -z (T c+k,.(tX.'(t' 
Gz2,s{k-t,t') = -i(t Ck,.(t)ag!;(t')) , wi^.t,3(k;t,t') = -i(t a%{t)c+,,{t') 



(3) 

where g = ia^ . Using the equations of motion for operators 0^0- and Ck.o- and transforming 

to the frequency representation, we obtain the following set of equations for the Green's 
functions: 

(a; - £,2(k)) G,2(k,a;) - iy3(k)G'+ 3(k, ^) + A,2F+(k,a;) = 1, (8) 

{lo - £3(k)) G'+_3(k,a;) - W^{k)GA^.^) + A3(k)F+_3(k, a;) = 0, (9) 

{u + e,2{k)) F+(k,cu) + 1^3(k)Ft^3(k,a;) + A%Gz2{k,uo) = 0, (10) 

(a; + £3(k))F+ (k,.;) + iy3(k)F+(k,u;) + A+(k)G'+3 (k,.;)=0, (11) 



where 



A3(k) = A33(k) = A,^ ^fi , n.r .2 ■ (12) 



{Bp - £3)2 + {W,r + iWy^' 
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In the following, we will make use of quasiparticle density of states of the band, 

N{uj) = --jlmGUuj,k)j^^. (13) 

The retarded Green's function is obtained by the solution of the system of equations (8-11) and 
has the form: 

'~ uj- E^(k) + iO^ uj- E+{k) +iO^ UJ + E_{k) + iO^ u + E+{k) + lO ' 
where the coherence factors are defined as 



,2 



{E^ + e.2){-El + ej + Aj) + WjjE. - £3) 
2E_{El-E^_) 

(i?_ - e,.){-El + el + Ai) + Wi{E. + e,) 



"'^''^ ^ 2E_{El-E^_) 
The dispersion relations for two branches of the excitation spectrum have the form 



± 



4(k) + - em - Ai(k))2 + 4l^|(k){(A.. - A3(k))2 + (£,.(k) + 8^^} 



The set of equations (8-11) can be reduced to the system of two equations for Green's 
functions 6*^2 and with the effective order parameter defined by the formula 

A3(k)iy|(k) , 

^^^^^) = £i(k) + Ai(k)+^^^- ^^^^ 

The sign of the order parameter A3(k)(12) is determined by the sign of the parameter Ad{7), 
which can be positive or negative depending on the type of interaction Ud (repulsion or attrac- 
tion). The parameter A^ = if Ud = or Ud > U^"^ > 0, where Uf^^ is a certain critical 
value of repulsion on the (1^2 -y2 orbital of copper which suppresses superconductivity in the 
system. As seen from the formula (15), in the case of a nonzero interaction Ud on the dx2_y2 
orbital (A^ 7^ 0), the order parameter depends on quasimomentum, and in the case of repulsion, 
< {Ud > 0), the parameter A22(k) changes its sign. 
The approach considered explains the strong anisotropy of the order parameter as a con- 
sequence of the symmetry properties of the one-particle matrix element of the interband hy- 
bridization PV3(k) which has nodes along the diagonals of the Brillouin zone due to the 
difference between the types of symmetry of the initial bands. In this case, the branch E_[k) 
of the excitation spectrum vanishes at points located on lines in the k space ^ at which the 
effective order parameter Az2{k) is equal to zero. It should be noted that the initial interaction 
is regarded as isotropic; i.e., the effect considered does not depend on the origin of the pairing 
mechanism. In addition, the approach does not require a strong anisotropy of the spectrum 
and, hence, can be used for various types of high-temperature superconductors. The model 
proposed in [T7j also leads to formulas of the type (14), (15), but it is based on the exotic 
condition of the interaction sign reversal in various regions of the Fermi surface. 
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The diagonalization of the one-particle part of Hamiltonian (1) leads to the problem with 
anisotropic attraction in energy bands. The effective order parameters in the bands may have 
nodes, but have no pure d or (s + d) symmetry since their anisotropy is determined by the band 
representations of the space symmetry group of the lattice [121 HH UH] ■ 

The formulated model of the order parameter anisotropy makes it possible to obtain the - — 

dV 

characteristics of SIS junctions for superconductors of the BSCCO type. Taking into account 
the crystal structure of these compounds, we will assume that tunneling along the c axis occurs 
mainly through the (ij,2 orbitals of copper in the Cu02 plane (and through the apical oxygen 
which is not included explicitly in the model). The matrix element of tunneling is assumed 
to be independent of momentum (Tkp = T = const) due to the random formation of bonds 
between the dz2 orbitals on both sides of the break junction. Thus, tunneling between two 
superconductors {Cu02 layers) is treated as occuring through a number of point junctions. 

The expression for the dependence of the quasiparticle tunnel current on the voltage applied 
to the junction in this case assumes the standard form: 

oo 

I{y) = 4e|Tp / [n{uj) - n{uj - eV)]N{uj)N{u - eV)du , (16) 



where n{u) = [eT +1]~^ is the Fermi distribution function and N{uj) is defined by formula (13). 
All calculations were performed for temperature T = 0. An important aspect of this approach 
is that N{u) is not the average density of states in the conduction band (as in the single-band 
BCS model), but is a partial density on the d^^ orbital. 

In order to find the density of states N{uj) and the tunnel conductance in accordance 

with the chosen approximation, we used the values of parameters from the tight-binding model, 
which correspond to the calculations of the band structure of HTSC [20 . If all parameters are 
expressed in terms of A^i and the position of energy levels is measured from the chemical 
potential, we have td-p = 75, t^i^p = 15, Sd = —75 and Sp = —50. In all calculations, we 
assumed the existence of a finite relaxation constant 7 = 0.05. The parameter determining the 
width of the initial £22(k) band is = 1.5. 

For the chosen values of the parameters, we analyzed the behavior of the characteristics N{uj) 

and depending on the interaction value between electrons on the rfx2-y2 orbital of copper 

and on the position of the center £^,2 of the initial z"^ band relative to the chemical potential 
(Figs. 2-4). The obtained dependencies were compared with these of single-band model with 
the s- and rf-types of symmetry of the order parameter (A = const and A ~ {cos{kx) — cos{ky)) , 
respectively). 

Figure 2A shows the A^(c<;) curves calculated for the normal state (A^2 = 0). These curves 
display split peaks associated with the Van Hove singularities of the initial 6x2 (k) band for 
various positions of its center relative to the chemical potential, as well as the following peaks 
formed by singularities of the £3(k) band. Figure 2b shows the same curves in the case when 
the £22 (k) spectrum is replaced by a dispersionless level (^^2 = 0). 

Fig. 3. A presents the N{u) dependencies in the superconducting state for the case when the 
2^ band is far from the chemical potential {e^^ = —10). It can be seen from the figure that if 
there is electron-electron repulsion on the d^'^-y'^ orbital A^ = —0.5), the density of quasipar- 
ticles is similar to that calculated in a model with the d symmetry of the order parameter (in 
particular, it is a linear function of uj in the vicinity of u; = 0). If there is no interaction on 




co/A m/A 

Figure 2: Normal density of states in the band: A) = 1.5, B) = at = (curves 
a), —1.5 (curves b) and —10 (curves c). Peaks 1 and 1' are due to the Van Hove singularities of 
the initial band, peaks 2 and 2' are formed by the Van Hove singularities in the £3(k) band, 
and peaks 3 are associated with the level. 

the dx2-y2 orbital, the distance between the superconducting peaks and their height increase. 
According to analytical estimations, N(io) ~ u^^"^ for very small frequencies. Similar results 
can also be obtained in the single-band model with a nontrivial spectrum and an effective order 
parameter of certain symmetry. As noted above, taking into account the above approximations, 
the diagonalization of the single-particle part of Hamiltonian (1) leads to the single-band model 
with anisotropic pairing. In this case, the symmetry of the corresponding order parameter can 
be approximated by an (s + d)-tjpe symmetry. 

If the center of the initial band is close to the chemical potential (see Figs. 3B and 3C), 
the difference between the a and b curves (with and without taking into account of repulsion) 
vanishes due to the dominant role of the Van Hove singularities on energy scales of the order 
of Az2. In both cases, the behavior of the density of quasiparticles becomes of the s type. The 
calculated density of states demonstrates the asymmetry of peaks, which is associated with the 
position of Van Hove singularities in the normal density of states (see Fig. 2). The Van Hove 
singularity nearest to the chemical potential increases the height of the corresponding peak in 
the density of states in the superconducting state. A comparison of Figs. 3B and 3C shows 
that, for a certain intermediate value of 6^2, there is the mirror switching of the asymmetry of 
peaks associated with the model spectrum structure. 

The curves presented in Fig. 3A show that even without the repulsion on the dj.2_y2- 
orbital if the ez2 level is far from the chemical potential {ez2 = —10), the distance between 
the superconducting peaks is smaller than 2A^2. As the ez2 level approaches the chemical 
potential, the distance between the superconducting peaks increases, and the position of the 
peaks on the N{uj) curve for 6^2 = corresponds to the parameter A^2. In this case, the 
height of the peaks increases significantly. The distance between the peaks changes as a result 
of the displacement of the ez2 level relative to the chemical potential for the fixed values of 
the remaining parameters. However, a more detailed analysis of this phenomenon requires the 
solution of the self-consistent equations for order parameters as functions of the parameters 
of the model as well as the determination of the interaction constants Uz2 and Ud taking into 
account the specific mechanism of pairing, which is beyond the scope of the present paper. The 
dependence of the distance between the superconducting peaks on the position of the chemical 
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co/A co/A 

Figure 3: Density of states in the normal and superconducting states in the band: A) 
£^2 = —10, t^i = 1.5, B) £^2 = —1.5, = 1.5, C) £^2 = 0, = 1.5, D) £^2 = —1.5, t^z = 0. 
Curves a are plotted in the absence of interaction on the d^^-yi orbital {Ud = 0); curves b take 
into account the repulsion on the d^z-yi orbital {Ud > 0); curves c plotted in a model with the 
s symmetry of the order parameter, curves d are plotted in a model with the d symmetry of 
the order parameter; curve e illustrates the normal density of states. 



potential in the model with the [s + (i)-type of order parameter symmetry was also found 
theoretically in [TB] . 

A comparison of the A^(a;) dependencies found in the model with single-band models with 
the s- and d-tjpes of the order parameter symmetry was carried out taking into account the 
location of the Van Hove singularities associated with the initial £22(k) band. This is pointed 
out by similar asymptotic behavior of densities of states calculated for both the normal and 
the superconducting state in s-, d- and considered models (see Fig. 3C) on the scales of energy 
higher than A22. 

A comparison of the results presented in Figs. 3B and 3D shows that, as the width of the 
initial z"^ band determined by the parameter tz^ decreases [t^^ = 1.5 in Fig.3B and t^2 = in 
Fig.3D), the heights of the superconducting peaks increases, since the singularity in the normal 
density of states in this case is displaced towards the chemical potential (see Fig. 2). 

Figure 4 shows the characteristics of SIS junction calculated by formula (16) using 
the determined N{uj) dependencies for three positions of the £^2 level relative to the chemical 
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Figure 4: characteristics for A) e 



-- -10, t^2 = 1.5, B) e,2 = -1.5, t,2 = 1.5, C) e,2 = 0, 

= 1.5, D) £^2 = —1.5, tz2 = 0. Curves a are plotted in the absence of interaction on the 
dx2-y2 {{Ud = 0); curves b take into account the repulsion on the dx2-y2 orbital {Ud > 0); curves 
c plotted in a model with the s symmetry of the order parameter, curves d are plotted in a 
model with the d symmetry of the order parameter. 



potential e^^ = —10, = —1.5, e^^ = 0, which can correspond to different HTSC types. The 
common feature of these characteristics is that in the case of repulsion the superconducting 
peak appears at lower voltages and is lowered, but the curve at low voltages lies higher than 
in the absence of interaction on the dx2-y2 orbital of copper. If the center of the initial 

band is far from the chemical potential (£^2 = —10), the characteristics are similar to those 

calculated in a model with the d symmetry of the order parameter (Fig. 4a). As the £^2 level 
approaches the chemical potential, the curves tend to the dependencies obtained in the model 
with the s symmetry of the order parameter (Figs. 4b and 4c). For low voltages V < A-^2, 
the IV characteristic of an SIS junction is strongly suppressed. In analogy with the A^(a;) 

dependencies, the - — characteristics exhibit extremely narrow superconducting peaks. This 
dV 

corresponds to a number of experimentally observed results for BSCCO-type compounds, e.g., 
optimally doped j^, overdoped 0, and underdoped [21] Bi2Sr2CaCu20s+s- 

Thus, the electron density of states and current-voltage characteristics of SIS junctions are 
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calculated in the framework of the multiband HTSC model with an anisotropic effective order 
parameter. We took into account the hybridization between the dx and dy orbitals of oxygen 
and the dx2_y2 orbitals of copper in the Cu02 plane and assumed the presence of electron- 
electron attraction on the (^^2 orbital of copper. The performed calculations demonstrate the 
possibility of explanation the results of experimental measurements of the electron density of 
states and the tunnel properties of SIN and SIS junctions (e.g., s-type characteristics of 
BSCCO-type compounds) on the basis of the model taking into account the structure of the 
electron spectrum of HTSC. It should be stressed that anisotropy of the order parameter and 
of the excitation spectrum is explained only by the symmetry of the crystal lattice and atomic 
orbitals and does not depend on the nature of pairing. 

The main result of this work is that for a strongly anisotropic order parameter A(k) (15) 
(e.g. of d-type) and in the presence of nodes in the excitation spectrum of a superconductor, the 
simple inclusion of the real band structure of the HTSC in a wide range of the model parameter 
values results in the experimentally observed [U El 121] s-type behavior of the current- volt age 
characteristics SIN and SIS junctions. A comparison of the curves calculated for various 
values of the parameters of the model and in the simplest cases of models with the s and d 
symmetries of order parameters indicates a strong dependence of both the density of states and 

the — characteristics on the electron band structure of high-temperature superconductors. 
dV 

These characteristics also depend on the crystal lattice symmetry and on the presence of an 
additional interaction between electrons, which was introduced in the model on the d^^^yi 
orbital of copper. Above analysis points to the necessity of considering of the band structure 
of HTSC for interpretation of the results of tunnel experiments . 

The authors are grateful to B.A. Volkov, E.G. Maksimov, and Ya.G. Ponomarev for nu- 
merous discussions in the course of the research. This work was supported by the Russian 
Foundation for Basic Research (N 02-02-16925, N 00-15-15998) and Landau Scholarship Grant. 
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